Probabilistic implementation of universal quantum processors 
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We present a probabilistic quantum processor for qudits. The processor itself is represented by a 
fixed array of gates. The input of the processor consists of two registers. In the program register the 
set of instructions (program) is encoded. This program is applied to the data register. The processor 
can perform any operation on a single qudit of the dimension A'' with a certain probability. If the 
operation is unitary, the probability is in general 1 /N'^ , but for more restricted sets of operators the 
probability can be higher. In fact, this probability can be independent of the dimension of the qudit 
Hilbert space of the qudit under some conditions. 



I. INTRODUCTION 

Schematically we can represent a classical computer 
as a device with a processor, which is a fixed piece of 
hardware, that performs operations on a data register 
according to a program encoded initially in the program 
register. The action of the processor is fully determined 
by the program. The processor is universal if we can real- 
ize any operation on the data by entering the appropriate 
program into the program register. 

In this paper we shall examine a quantum version of 
this picture. Specifically, in close analogy with recent 
papers by Nielsen and Chuang [|| and Vidal and Cirac 
1^, we will study how a quantum program initially put 
into a program register can cause a particular operation 
to be applied to a data register initially prepared in an 
unknown state. We shall first consider the case in which 
the data consists of a single qubit, and the program of 
two qubits. We shall then examine higher-dimensional 
systems. 

Nielsen and Chuang [Q originally formulated the prob- 
lem in terms of a programmable array of quantum gates, 
which can be described as a fixed unitary operator, P^p, 
that acts on both the program and the data. The ini- 
tial state, \^u)p, of the program register stores informa- 
tion about the one-qubit unitary transformation U that 
is going to be performed on a single-qubit data register 
initially prepared in a state The total dynamics of 

the programmable quantum gate array is then given by 



(1.1) 



where only pure data states were considered. The pro- 
gram register at the output of the gate is in the state 
|2[/)p - which was shown to be independent of the input 
data state 

Nielsen and Chuang proved that any two inequivalent 
operations U and V require orthogonal program states, 
i.e. — 0. Thus, in order to perfectly imple- 

ment a set of inequivalent operations, {[/_,■ |j G J}, the 



state space for the program register must contain the or- 
thonormal set of program states, {|S(7.)|j G J}. This 
means that the dimension of the program register must 
be at least as great as the number of unitary operators 
that we want to perform. Since the set of unitary opera- 
tions is infinite, the result of Nielsen and Chuang implies 
that no universal gate array can be constructed using fi- 
nite resources, that is, with a finite dimensional program 
register. They did show, however, that if the gate ar- 
ray is probabilistic, a universal gate array is possible. A 
probabilistic array is one that requires a measurement to 
be made at the output of the program register, and the 
output of the data register is only accepted if a particu- 
lar result, or set of results, is obtained. This will happen 
with a probability, which is less than one. 

Vidal and Cirac |^ have recently presented a proba- 
bilistic programmable quantum gate array with a finite 
program register, which can realize a one parameter fam- 
ily of operations, where the parameter is continuous, with 
arbitrarily high probability. The higher the probability 
of success, the greater the dimensionality of the program 
register, but the number of transformations that can be 
realized is infinite. They have also considered approx- 
imate programmable quantum gate arrays, which per- 
form an operation Ejj very similar to the desired U , that 
is F{Eu, U) > 1 — e for some transformation fidelity F. 

Another aspect of the encoding of quantum operations 
in the states of program registers has been discussed by 
Huelga and coworkers |^] . In this paper the implementa- 
tion of an arbitrary unitary operation U upon a distant 
quantum system has been considered. This so called tele- 
portation of unitary operations has been formally repre- 
sented as a completely positive, linear, trace preserving 
map on the set of density operators of the program and 
data registers: 



(1.2) 



\'E.u)p(g>\lp)d] = \^u)ap(^{U\lp)d) 

Here |^)a represents a specific entangled state that is 
shared by two parties, Alice and Bob, who want to tele- 
port the unitary operation U from Alice to Bob. Huelga 
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et al. g have investigated protocols which achieve the 
teleportation of U using local operations, classical com- 
munication and shared entanglement. 

In the present paper we will address the problem of 
implementing an operation U , encoded in the state of a 
program register |Sc/)p, on the data state 1-0)^- The gate 
arrays we present are probabilistic; the program register 
must be measured at the end of the procedure. In Sec- 
tion H we present a simple example of how to apply an 
arbitrary operation to a single qubit initially prepared in 
a state The gate array consists of four ControUed- 
NOT (C-NOT) gates, and can implement four programs 
perfectly. These programs cause the one of the opera- 
tions 1, <Tx, —icTy, or ffz to be performed on the data 
qubit. Here 1 is the identity and Uj, where j = x,y,z 
is a Pauli matrix. By choosing programs that are lin- 
ear combinations of the four basic ones, it is possible 
to probabilistically perform any linear operation on the 
data qubit. In Section [II we generalize the idea to an 



arbitrary dimensional quantum system, a qudit. 



II. OPERATIONS ON QUBITS 

We would like to construct a device that will do the 
following: The input consists of a qubit, \^)d, and a sec- 
ond state, |S(7)p, which may be a multiqubit state, that 
acts as a program. The output of the device will be a 
state U\'ip)d, where U is an operation that is specified by 
In order to make this a little less abstract, we first 
consider an example: Let 10) and \4>±) be two orthogonal 
qubit states, and suppose that we want to perform the 
operation 



>(0x|-|0)(0| =1-2|0)(0|, 



(2.1) 



on \ip)d- The action of this operator is analogous to that 
of (Tz in the basis {|0), |1)}, except that it acts in the basis 
110^), 10)}. That is, (Tz does nothing to |0) and multiplies 
|1) by —1, while does nothing to |0^) and multiplies 
10) by —1. Can we find a network and a program vector 
to implement this operation on |V')d? 

We can, in fact, do this by using the network for a 
quantum information distributor (QID) as introduced in 
Ref. [Q (this is a modification of the quantum cloning 
transformation [^,^ ). In this network the program reg- 
ister is represented by a two qubit state |5yi)p. Before 
we present the network for the programmable gate ar- 
ray, we shall introduce notation for its components. A 
ControUed-NOT gate Djk acting on qubits j and k per- 
forms the transformation, 



Djk\m)j\n)k = \m)j\miSn)k, 



(2.2) 



where j is the control bit, k is the target bit, and m and 
n are either or 1. The addition is modulo 2. The QID 
network consists of four ControUed-NOT gates, and acts 
on three qubits (a single data qubit denoted by a sub- 
script 1 and two program qubits denoted by subscripts 2 



and 3, respectively). Its action is given by the operator 
^123 — D31D21D13D12. As our first task, we shall deter- 
mine how this network acts on input states where qubit 
1 is in the state jV'), and qubits 2 and 3 are in Bell basis 
states. The Bell basis states are defined by 
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We find that 







^'l23lV')l|*+)23 = (a,|^)i)|$ + ); 

A23lV')l|^ + )23 = l^)l|*+); 
Pl23|V')l|*-)23 = (az|V)l)|«'-). 



(2.3) 



(2.4) 



Any operation on qubits can be expanded in terms of 
Pauli matrixes and the identity. The above equations 
mean that the Bell basis vectors are "programs" for a 
complete set of operations. In order to see how to make 
use of this, let us expand our proposed operation in terms 
of this complete set. Expressing |0) as |0) = /i|0) + i/|l), 
we have that 



A, = 1-210) (01 = 



= — (/ii^* + fj,*iy)<7x + (fiiy* — ii*v){—iay) 



(2.5) 



We can now apply the operation A to 
the "program" vector 



/iV)|$H 



/23 
)23, 



by sending in 
(2.6) 



and measuring the program outputs in order to deter- 
mine if they are in the state (|$+) + |$_) -I- \^^))/V3. 
If they are, our operation has been accomplished. Note 
that the measurement is independent of the vector |0) 
so that no knowledge of this vector is necessary to make 
the measurement and to determine whether the proce- 
dure has been successful. As we see, the probability of 
success is 1/3 for the implementation of the operation 
Az which is parameterized in general by two continuous 
parameters (i.e. the state |0)). 

Let us examine the program vector more carefully. If 
we define the unitary operation, Uinit, by 



C/™t|00) = -|10) 
u^mt\n) = |01); 



C/^mdlO) =-|ll) 

U^mtlOl) = |00); 



(2.7) 
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we have that 



so that 



|Sa)i2 = U„,u^mcj,^) + (2.8) 

Finally, we can summarize our procedure. The steps 
are 

1. Start with the state -i=(|0)|0i) + \(l)±)\(f>)) . 

2. Apply Uinit- 

3. Send the resulting state into the control ports (inputs 

2 and 3) and 1-0) into port 1. 

4. Measure (|$+) + |<I>_) + |*_))/\/3 at the output of 

the control ports. 

5. If the result is yes, then the output of port 1 is 

Before proceeding to a more general consideration of 
this network, let us make an observation. Suppose that 
we carry out the same procedure, but instead of starting 
with the program vector {\(f>)\(f)±) + \(f)±)\(f))) /^/2, we start 
instead with the program vector (|0)|0) — \(j)±)\4'±))/V2. 
At the end of the procedure the output of the data reg- 
ister is 1-0), where 



1- E i^^'^i' 



(2.13) 



Now let us go back to our network and consider the 
program vector given by 



j,k=0 



(2.14) 



and at the output of the program register we shall mea- 
sure the projection operator corresponding to the vector 
(1/2) J2j k=o I'^jk)- If the measurement is successful, the 
state of the data register is, up to normalization, given 

by 



(2.15) 



After this state is normalized, it is just (1/|| A0|j)|'!/')- 
This means that for any transformation of the type given 
in Eq. ( p. 10 ) , we can find a program for our network that 
will carry it out. 




(2.9) 



The operation interchanges |0) and \4>±)- Its action 
is analogous to that of ax , which interchanges the vectors 
|0) and |1). The probability of success for this procedure 
is also 1/3. 

We now need to determine whether there is a program 
for any operator that could act on 1-0). The operator need 
not be unitary; it could be a result of coupling ji/i) to an 
ancilla, evolving the coupled system (a unitary process) , 
and then measuring the ancilla. Therefore, if A is now 
any linear operator acting on a two dimensional quantum 
system, the transformations in which we arc interested 
are given by 



\A^\ 



(2.10) 



Let us denote the operators, which can be implemented 
by Bell state programs, by 6*00 = 1, <5'oi = o^x, •S'lo = ^zj 
and Sii = —lay. Any 2x2 matrix can be expanded in 
terms of these operators, so that we have 



j,k=0 

We now define ajk — (ijk/ where 



j,k=0 



(2.11) 



(2.12) 



III. GENERALIZATION TO QUDITS 

In order to extend the network presented in the previ- 
ous section to higher dimensions, we must first introduce 
a generalization of the two-qubit C-NOT gate jj] (see also 
Ref. §). As we noted previously, it is possible express 
the action of a C-NOT gate as a two-qubit operator of 
the form 



D„ 



k,m—0 



\k)a{k\ (g) \m(Bk)b{m\ 



(3.1) 



In principle one can also introduce an operator de- 
fined as 



Dib= J2 \k)a{k\®\mek)b{m\ 

k,m—0 



(3.2) 



In the case of qubits these two operators are equal, but 
this will not be the case when we generalize the opera- 
tor to Hilbert spaces whose dimension is larger than 2 
In particular, we can generalize the operator D for 
dimension A'^ > 2 by defining 

Dab^ J2 \k)a{k\<»\{m + k)modN)h{m\ , (3.3) 

k,m—Q 

which implies that 
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N-1 



|(m — fc)modiV)f,(r7i| 



(3.4) 



k,m—0 



From this definition it follows that the operator Dab acts 
on the basis vectors as 



Dab\k)\m) = |/c)|(fc + m)modiV) 



(3.5) 



which means that this operator has the same action as the 
conditional adder and can be performed with the help of 
the simple quantum network discussed in . Now we see 
that for iV > 2 the two operators D and do differ; they 
describe conditional shifts in opposite directions. There- 
fore the generalizations of the C-NOT operator to higher 
dimensions are just conditional shifts. 

In analogy with the quantum computational network 
discussed in the previous section, we assume the network 
for the probabilistic universal quantum processor to be 



^123 — 



(3.6) 



The data register consists of system 1 and the program 
register of systems 2 and 3. The state \^u)23 acts as the 
"software" which the operation to be implemented on the 
qudit data state l^*)!. The output state of the three qu- 
dit system, after the four controlled shifts are applied, 
reads 



\Q) 



123 



^3lC|l£'l3£'l2|*)l|Sa)23 



(3.7) 



A graphical representation of the logical network (|3.7| ) 
with the conditional shift gates Dab in Fig- 0- 
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FIG. 1. A logic network for the universal qua ntum proces- 
sor as given by the unitary transformation (3.7). The action 
of the controlled shift operator Djk is represented as follows: 
The control qudit is represented by • while the target qu- 
dit is represented ffi with the rightarrow. The action of the 
operator Dj^ is represented by left arrow. 

The sequence of four operators acting on the basis vec- 
tors gives |ri)i|m)2|fc)3 as 



D3iDl^Di3Di2\ n)i|rn)2|fc)3 = \{n — m + k)modN)i \{m + n)mod N)2 \{k + n)mod A^)j 



(3.8) 
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We now turn to the fundamental program states. A ba- 
sis consisting of maximally entangled two-particle states 
(the analogue of the Bell basis for spin-i particles) is 
given by [§ 

|S„j„) = ^= E exp(^i-^m/c^|A:)|(fc - n)modiV) , (3.9) 

where m, n = 0, . . . , TV — 1. If |Sm„)p is the initial state 
of the program register, and l^f) = X^j'^ib)^ (here, as 
usual, J2j lo^iP = 1) is the initial state of the data regis- 
ter, it then follows that 



-pL23|^')l|>^mri)23 — 

Ea-j 2'Kikm i.,,,,,, , 

-j=^^V^^Pi^sb)\k)\k-n) 



N 
2Trikm 



N 



\j - n)\k + j)\k + j - n) 



jk 

/ -1= exp ■ 

E—2iTijm. . M„ , 

Jk 

([/("")|vE'))|S„ 



TV 



(3.10) 



where we have introduced the notation 



N-l 

E 

s=0 



exp ■ 



-2iTTsm . 



n){s\ 



(3.11) 



This result is similar to the one we found in the case 
of a single qubit. We would now like to examine which 
transformations we can perform on the state in the data 
register by using a program consisting of a linear com- 
bination of the vectors |Sm„) followed by the action of 
the processor P123 and a subsequent measurement of the 
program register. 

The operators [/("'") satisfy the orthogonality relation 



Tr 



([/(™'n'))t[/(n"0 



(3.12) 



The space of linear operators T{TC) defined on some 
Hilbert space Ti. with the scalar product given by (3.12) 
we know as Hilbert- Schmidt space. Thus the unitary op- 
erators f/^™") form an orthogonal basis in it and any 
operator A e T{T-l) can be expressed in terms of them 



A 



N-l 

E 9" 

m,n— 



(3.13) 



The orthogonality relation allows us to find the expansion 
coefficients in terms of the operators 
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N 



Tr 



A 



Equations ( [3.12 ) and (3.13) imply that 



E 



= ±Tr{A^A). 



N 



(3.14) 



(3.15) 



Therefore, the program vector that implements the op- 
erator A is given by 



\va)23 



N 



Tr(AtA) 



1/2 N-1 

'Zmn|Sm„)23- (3.16) 

?n,n— 



Application of the processor to the input state |^')i|wa)23 
yields the output state 

I^^)l23 = E ^-nf^^""^!*)! ® |2™n)23- (3-17) 



Af-1 



|0*)=^/3,*|A;). 



(3.23) 



fc=0 



We then have that 



|$)23 = (W^2®l3)(i?L)'(|2 



00/23 



ir)2i</')3 



(3.24) 

A network that performs the operation (1^2^13) (j^23^ 
could be added to the input of the program register so 
that the si mpler state that appears on the right-hand 
side of Eq. ( 3.24 ) could be used as the program. At the 
output of the processor we have to perform the projec- 
tive measurement discussed in the previous paragraph, 
and the probability of achieving the desired result is the 
same as the probability of successfully implementing the 
transformation, A. In this case the probability is l/N"^. 



To obtain the final result we perform a projective mea- 
surement of the program register onto vector |M)23 



(3.18) 



m.n— 



If the outcome of the measurement is positive, then we 
get the required transformation A acting on an unknown, 
arbitrary input state l^")!. 

Let us consider an example. Suppose we choose for A 
the unitary operator 1 — 2|(/))((/)|, where the normalized 
state \4>) can be expressed as 



N~l 



(3.19) 



fc=0 



The expansion coefficients for this operation are given by 



N-l 



<5mO<^nO 



N 



J2 e2"^"/^/3*/3fc_„, (3.20) 



fc=0 



and the program vector for this operation is 

N-l 



1^)23 — 1^00)23 



-= J2 filkP-ik+n)\k)2\k- 
k,n—0 



n)3- 



(3.21) 



The program vector can be obtained from a state more 
closely related to \(f)) if we introduce a new unitary oper- 
ator and a "complex conjugate" vector. Define the oper- 
ator W by 



W\k) = \-k), 
and the vector by 



(3.22) 



IV. SUCCESS PROBABILITY 

The probability, p, of successfully applying the opera- 
tor A to the state in our example is rather small. 
This is because the operator we chose was a linear com- 
bination of all of the operators This means that 
if the data register consists of I qubits, i.e. TV = 2', then 
the probability of a successful implementation of a gen- 
eral transformation A decreases exponentially with the 
size of the data register. However, if we were to choose 
an operator, or set of operators, that was a linear com- 
bination of only a few of the then the success 
probability can be significantly improved. This would en- 
tail making a different measurement at the output of the 
program register. Instead of making a projective mea- 
surement onto the vector |M), one would instead make 
a measurement onto the vector 



J^l/2 E I 

7n,n:qmn7^0 



(4.1) 



where Af is the total number of n onzero coefficients qmm 
in the decomposition in Eq. ( [3.13| ) . If the operation being 
implemented is unitary, then, in this case, the probability 
of implementing it is 



P 



1 



(4.2) 



where M is the total mmrber of nonzero coefficients qmm 
in the decomposition (3.13). There are, in fact, large 
classes of operations that can be expressed in terms of 
a small number of operators C/^™"^ ||l^. For these op- 
erators, the probability of success can be relatively large 
and, in principle, independent of the size of the Hilbert 
space of the data register. 
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Example 1. 

Let us consider the one-parameter set of unitary trans- 
formations 



U^p = cos ipt + i sin ip 



1 + i 



[/(Ol) 



1 - i 



[/(03) 



(4.3) 



where the unitaries C/^™"-' are given by Eq.(3.11). These 
unitaries for = 4 can be exphcitly written as 

[/(oi) ^ ; C/(03) ^ , (4.4) 

where Pg — \ s) (s \. From here we find the expression for 
the operator ( |4.3| ) in the form: 

= cosifit + i sin ip [Pq + Pi - P2 - P3] , (4.5) 

We note that if we rewrite the parameters s as binary 
numbers, s — ji2 + jo, where jk is either or 1, and 
express the states \s) as tensor products of qubits, i.e. 
|s) = bo) J we find that the operator in brackets on 

the right-hand side of Eq. (E^) can be expressed as 



1 + i 



C/(oi) 



1 - i 



[/(03) 



era «) 1. 



(4.6) 



From Eq. ( |4.5[) it is clear that U^p has eigenvalues of mag- 
nitude 1, which implies that U^p is unitary. It can be 



realized by the universal quantum processor (3^) with 
a probability of successful implementation equal to 1/3. 
This example illustrates that it is possible to realize large 
classes of unitary operations with a probability that is 
greater than the reciprocal of the dimension of the pro- 
gram register. 

This example can be easily generalized. Consider 
a one-parameter set of unitary operators acting on a 
Hilbert space consisting of I qubits, which is given by 



Uip = cos (pi + i sin p 0-3 



1 



(4.7) 



The operator 0-3 (g) is diagonal and therefore only 

the diagonal unitaries from our s et C/^™"-', i.e. C/(™o), 
appear in its expansion, Eq. (3.13) . Moreover the coef- 
ficients Qmo in the expansion arc non-vanishing only for 
odd TO. It follows that 



C/^ — cos(p 1*^ 



odd rn 



.0) 



(4.8) 



and the probability of a successful implementation of this 
unitary transformation is p 2/(2' + 2). 
Example 2. 

For some sets of operators it is possible to do even better 
than we were able to do in the previous example. Con- 
sider the one-parameter set of unitary operators given 
by 



Ui) = cosi31 + i sin d ;7(0'^/2) , 



(4.9) 



where N is assumed to be even. That this operator is 
unitary follows from the fact that f/(0'^/2) is self-adjoint. 
A program vector that would implement this operator is 

|$)23 = COSl?|Soo)23 -I- isini?|So,Ar/2)23, (4-10) 

and at the output of the program register we make a 
projective measurement corresponding to the vector 

1^)23 = -^(|Soo)23 + |So,Ar/2)23- (4.11) 

The probability for successfully achieving the desired re- 
sult, i.e. the vector U^\'^)i in the data register, is 1/2 
irrespective of the value N, i.e. the number of qubits. 



V. CONCLUSION 

We have presented here a programmable quantum 
processor that exactly implements a set of operators 
that form a basis for the space of operators on qudits. 
This processor has a particularly simple representation 
in terms of elementary quantum gates. It is, however, 
by no means unique. It is possible, in principle, to build 
a processor that exactly implements any set of unitary 
operators that form a basis for the set of operators on 
qudits of dimension N, and uses any orthonormal set 
of N"^ vectors as programs. Explicitly, if the set of op- 
erators is {Ki|n — 1,...N^} and the program vectors 
are {|yn)|n = 1, . . . N^}, the processor transformation is 
given by 



Pdn = 



E 



(5.1) 



where the superscript (d) on the operator Vn indicates 
that it acts on the data register. 

As an example, consider a data register consisting of 
I qubits. We could use the processor discussed in sec- 
tion III to perform operations on states in this register, 
but we can also do something else; we can use I single- 
qubit processors, one for each qubit of the data register. 
Specifically, our unitary basis for the set operations on 
the data register would be 



UjK = Uj^ku...:3iki 



(5.2) 



where J = {ji, ■ ■ ■ ,ji) and K — {ki . . . ,ki) are sequences 
of zeros and ones, and the operators Sj^k^ are the de- 
fined immediately after Eq. (2.10). The program register 
would consist of I pairs of qubits, 21 qubits in all, with 
each pair controlling the operation on one of the qubits 
in the data register. Each of the operators in our basis 
can be implemented perfectly by a program consisting of 
the tensor product state, nL=o 1"^™!:™)' '^'^sre [S^™!,^) 
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is a two-qubit state that implements the operation Sj^^k,„ 
on the mth qubit of the data register. 

We are then faced with the problem of which proces- 
sor to use. This very much depends on the set of op- 
erations we want to apply to the data. How to choose 
the processor so that a given set of operations can be im- 
plemented with the greatest probability, for a fixed size 
of the program register is an open problem. A second 
issue is simplicity. One would like the processor itself 
and the program states it uses to be as simple as pos- 
sible. The simplicity of the processor is related to the 
number of quantum gates it takes to construct it. We 
would maintain that the processors we have presented 
here are simple, though whether there are simpler ones 
we do not know. Judging the simplicity of the program 
states is somewhat more difficult, but they should be re- 
lated in a relatively straightforward way to the operation 
that they encode. In many cases these states will have 
been produced by a previous part of a quantum algo- 
rithm, and complicated program states will mean more 
complexity for the algorithm that produces them. The 
program states proposed by Vidal and Cirac and the ones 
proposed by us in section II are, in our opinion, simple. 

A final open problem that we shall mention, is finding 
a systematic way of increasing the probability of success- 
fully carrying out a set of operations by increasing the 
dimensionality of the space of program vectors. Vidal 
and Cirac showed how to do this in a particular case, 
but more general constructions would be desirable 
Doing so would give one a method of designing programs 
for a quantum computer. 
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